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Abstract 

This work investigates an upper-limit of charge for a black hole in a nonlinear Newtonian theory of 
gravity. The charge is accumulated via protons fired isotropically at the black hole. This theoretical 
study of gravity (known as 'pseudo-Newtonian') is a forced merger of special relativity and Newtonian 
gravity. Whereas the source of Newton's gravity is purely mass, pseudo-Newtonian gravity includes 
effects of fields around the mass, giving a more complete picture of how gravity behaves. Interestingly, 
pseudo-Newtonian gravity predicts such relativistic phenomena as black holes and deviations from 
Kepler's laws, but of course, provides a less accurate picture than general relativity. Though less 
accurate, it offers an easier approach to understanding some results of general relativity, and merits 
interest due to its simplicity. The method of study applied here examines the predictions of pseudo- 
Newtonian gravity for a particle interacting with a highly charged black hole. A black hole with a 
suitable charge will reach an upper limit (expressed by pseudo-Newtonian gravity) in charge capacity 
before Coulomb's law repels like-charge particles away from the hole. In particular, this work attempts 
to push pseudo-Newtonian gravity to its extreme and discover how its results differ from general 
relativistic predictions involving the same proton bombardment. It is found that the results for an 
upper limit of charge in general relativity and this nonlinear theory of Newtonian gravity differ by 
a factor of four. This may give insight into the importance of space-time curvature effects on the 
description of particle dynamics around a black hole. 
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I. INTRODUCTION 



The fundamental law of gravitation discovered by Newton has withstood more than three 
centuries of experimental scrutiny and, at least in the weak gravity region throughout our 
Solar System, has provided an excellent description of gravitational phenomena. In terms of a 
field theory, Newtonian gravity is a linear theory completely described by a single scalar field 
variable, 0at, whose single source is mass. 

The comprehensive gravitational theory proposed by Einstein depicts gravitational phenom- 
ena as the result of a particular curvature of four- dimensional space. Any such curvature may 
result from a variety of causes, such as mass, fields, rotations of masses, etc. The complex 
mathematical nature of Einsteinian gravity as expressed by a system of coupled, nonlinear field 
equations makes obtaining any complete description of a specific curvature extremely difficult. 
The rewards of successful efforts, however, have led to some drastic differences from Newtonian 
gravity in cases of extremely strong gravitational fields. The most notable case, completely 
absent from Newtonian gravity, is that of gravity in the immediate vicinity of a black hole. 

It was shown by Peters^ that extending the basic notion of the source of Newtonian gravity 
to include fields as source contributions results in a nonlinear field theory, but one in terms of a 
single scalar field variable. Solutions in this "pseudo-Newtonian" gravity representing fields in 
regions exterior to a highly compacted mass, M, are found to possess singularities which have 
black hole characteristics. It was later shown by Young 2 that if M were imagined to have a 
net charge, Q, a black hole-like gravitational field would result in which the charge contributed 
directly and significantly to that gravitational field. Peters' solution bears a strong similarity 
to the well-known Schwarzschild solution found in Einsteinian gravity, but has the advantage 
of an enormous mathematical simplicity. Young's solution, on the other hand, contains certain 
features exhibited by the Reissner-Nordstrom solution of the Einstein-Maxwell equations, but 
again with the substantial mathematical simplicity offered by the pseudo-Newtonian formula- 
tion. 

The gravitational field of a charged black hole leads to new issues in the dynamics of charged 
particles interacting with the field. For example, a particle of mass m and charge q would ex- 
perience the competing effects of gravitational mass-attraction and Coulomb charge-repulsion. 
If we imagine starting with an electrically neutral black hole and gathering protons, the re- 
sulting charge build-up within the hole would have the effect of eventually repelling additional 
protons. If this were not so, the hole could ultimately become "over-charged" and consequently 
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destroyed, as discussed in the next sections. The work discussed here will be directed toward 
determining the maximum charge that proton accretion could add to a black hole. Section II 
reviews and verifies the nonlinear Newtonian theory, while sections III and IV find low energy 
and high energy proton bombardment charge limits, respectively. Section V compares the low 
energy case to the high energy case, while Section VI derives and compares the general rela- 
tivistic result with the nonlinear Newtonian answer. Motivated by the chance to complement 
the pedagogical nature of pseudo-Newtonian gravity, we have used S.I. units. 



II. PSEUDO-NEWTONIAN GRAVITY 

The Newtonian gravitational field due to localized mass distribution p m is completely de- 
scribed by the potential function, fa, satisfying the linear Poisson equation 

W 2 fa = AnGp m . (1) 

Exterior to spherically symmetric p m representing total mass M, this has solution 

fa(r) = . 2 

r 

This potential is independent of charge and angular momentum. The sole source of Newtonian 
gravity is mass. 

Peters^ examined the results of including the effective mass of the field produced by and 
surrounding the mass, by the following. There is an energy density stored in the gravitational 
field; this is given by-2 

Taken in conjunction with the mass-energy relationship, E = mc 2 , the energy density can be 
interpreted field-mass density, 

If the mass density, the source of Newtonian gravity, were taken to include that associated with 
the field (p m — > p m + pf), (OQ) and (j4j) would combine to yield the nonlinear field equation 

V 2 0+^(V0) 2 = 4vrGp m . (5) 



This has a spherically symmetric solution, exterior to the spherical region of M, 

GM\ 



r:.U) .-. :2c 2 In ( 1 - ^ (6) 



The Peters potential, Eq. ([6]), has some features worth discussion. First, careful examination 
shows, as is necessary with all relativistic results, that 0(r) — > 0jv( r ) in the limit of infinite c. 
Further, for large r, there is no distinction between <p(r) and 0jv(r). If a body of mass M were 
imagined to be extremely compact, having radius R < from its center, it would then be 
confined to an infinitely deep potential well. The pseudo-Newtonian gravity proposed by Peters 
thus leads, in a simple mathematical way, to the existence of a black hole which is strikingly 
analogous to the Schwarzschild black hole of Einsteinian gravity (but completely absent from 
Newtonian gravity). 

The all-inclusive character, in terms of sources, of Einsteinian gravity allows black holes other 
than the Schwarzschild type. The first such solution to the combined Einstein-Maxwell field 
equations representing a spherically symmetric mass containing a spherically symmetric charge 
was found independently by Reissner and Nordstrom. The pseudo-Newtonian counterpart of 
this solution^ is realized by including the electric field energy density, 

u e = |(V^) 2 , (7) 



where Eq is the permittivity of a vacuum and if is the electrostatic scalar potential produced 

Q 

4ireor ' 



by the charge, viz, . Q . The field equation in this case where p m has the value of zero, is 



V 2 + ^W) 2 -^(V^) 2 = O. (8) 

This can be shown to have a spherically symmetric solution, representing the gravitational field 
of a mass M carrying uniform charge Q, given by 



0(r; M, Q) = 2c 2 In 

where 



cosh ( b—^\ — a— sinh (b— 



(9) 



a^^G; IeI/A. ( 10 ) 

The potential given in Eq. ([9]) is only ostensibly dominated by the charge, and is insensitive 
to the sign of that charge. It has a limiting form identical to the Peters potential as Q — > 0. 
Further, it also implies the possible existence of a black hole since exhibits a singularity as 
the argument of the logarithm becomes zero. This is the case, for given Q and M, at r = Rbh, 
where 

Rbh = 2b ,3 0n . (11) 



In 
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In order that such a radius exist, the argument of the logarithm must be positive. This condition 
requires that such a charged black hole can exist in pseudo-Newtonian gravity only if the mass 
is not over-charged. It is convenient to express this condition with the inequality 



Q 
M 



2 



< Aire G ~ 10 



-20 



c 

kg 
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(12) 



Rather remarkably, this is precisely the same condition found in the Reissner-Nordstrom 
solution. 4 Should the condition be violated, no black hole exists. In the case of a proton, 
the charge to mass ratio has a square value approximately 10 16 [C /kg} 2 . Thus, as has been well 
documented, a proton is not a black hole! 

Graphical comparison of the charged and uncharged potentials is facilitated, as is the ensuing 
analysis, by first carrying out a rescaling. The unitless parameters A and a are defined as 

r 1 Q 

X ~Wh a ~7^GT M (13) 

where Rj,h = GM/2c 2 , i.e. the uncharged pseudo- Newtonian black hole radius (also the general 
relativistic spherical black hole radius measured in isotropic coordinates).- The inequality given 
in Eq. (Ti~2l restricts a values to the range < o < 1, the upper limit corresponding to a highly- 
charged mass. In these terms, the potentials of masses with and without charge can be written 
as, respectively, 



^<P(r;M,Q)= In 



cosh ( — j sinh ^— 



^<f>(r; M) = In (l - ± ) . (14) 



III. NON-RELATIVISTIC PARTICLE DYNAMICS 

A particle of mass m and charge q interacting with compacted mass M having charge Q will 
interact gravitationally and electrically with the latter. The particle's mass will interact with 
4>(r; M, Q) and its charge with the electrical potential produced by Q. The electrical potential 
here is not altered from its elementary form by including energy density effects. The inclusion 
adds only to total mass, altering the gravity because of the mass-energy relationship. Since an 
equivalent charge-energy relationship does not exist, the equations governing pseudo-Newtonian 
gravity remain coupled and asymmetric .- 

Any such particle would have non-relativistic (v < c/10) total energy E = T + V , where 

aO 

V(r;M,Q)=m(j)(r;M,Q) + -^-. (15) 

Aneor 
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In the course of its motion, the particle would undergo the competitive effects of gravitational 
attraction (since < 0) and, for qQ > 0, Coulomb repulsion. Since the fields in which the 
particle is considered here are central (i.e., functions of r only), its total energy will be conserved 
throughout its motion. In the case of the particle accelerated toward M, Q from infinity with 
kinetic energy E, its kinetic energy at any location will be 



In a plot where E is a constant and V(r; M, Q) is superimposed, the difference between these 
curves at any r would correspond to the particle's kinetic energy at that point. Physically 
realistic motion can occur only where the kinetic energy is positive, or at least zero; more 
specifically, it cannot be negative. Equation ([16]) indicates that the kinetic energy will decrease 
to zero at the critical distance, r c . This results in the equality 



In the case of a particle approaching M, Q from infinity, where V = 0, this equality would 
correspond to an intersection of the superimposed total and potential energy curves and indicate 
a turning point, or point of closest approach, in its motion. 

The main purpose of the present work is to examine the process of proton accretion by a black 
hole and, in particular, the charge buildup within the hole due to that process. An extremely 
large value of Q on M leads to no black hole at all. We imagine starting out with a black 
hole having little, or even no charge, then firing in protons isotropically to preserve spherical 
symmetry. The value of Q thus increases as the process continues and so then will the size of 
the Coulomb repulsion experienced by subsequent protons. Can protons continue to be added 
indefinitely, thus destroying the black hole, or will the Coulomb effect eventually dominate the 
gravitational attraction once some value, Q ma x, is reached? The addition of protons would 
have the additional effect of increasing the mass of the hole; that addition, however, would be 
negligible in comparison to the size of any significant M. This investigation has been restricted, 
for now, to the case of radially-approaching, non-relativistic protons having total energy less 
than about 4.5 MeV (v < c/10). 

The kinetic energy of an approaching proton would change with distance from the hole. 
This could increase or decrease depending on the amount of charge, presumably accumulated 
by earlier protons, already captured by the hole. Gravitational attraction dominates for small 
Q, allowing all protons with < E < 4.5 MeV to be captured. As Q increases, the Coulomb 
repulsion of subsequent protons grows in size and reaching Rbh becomes more difficult. 



T — E — V(r;M,Q). 



(16) 



E = V(r c ;M,Q). 



(17) 
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The analysis of the motion is most easily carried out by means of the general dynamical 
considerations discussed above. Writing Eq. ( Tl6|) explicitly, gives 

J- = - In (cosh (^) - 1 sinh (*)) - * ° (18) 

2mc 2 2mc 2 \ \\J a \\J J m^^G A V ' 

For protons, the last term on the right side of Eq. (18) is approximately 1.1 x 10 18 ^. The 
requirement that T > means that an approaching proton will encounter a turning point in 
its motion at the value of A that causes T to become zero. 

The kinetic energy, T, is just the difference between a horizontal (constant) E/2mc 2 curve 
and a V(X; M, Q)/2mc 2 curve for differing values of charge. The value of E is chosen as the 
approximate maximum for non-relativistic protons for which E/2mc 2 has a value 1/400. No 
intersection occurs in the case of small Q, meaning that any proton having this energy (or, in 
fact any E > 0) would get into the hole. As charge is added, the scaled-potential energy curves 
increase in size and an intersection can occur. The critical value of Q is that which maximizes a 
V curve at value 1/400. Any larger Q will then cause the V curve to intersect with the constant 
energy curve, indicating a turning point at the location of the intersection. This illustrates that 
a definite limit on Q exists such that no additional protons at this energy will be able to gain 
entry into the black hole. Hence, eventually Coulomb will become dominant. It is impossible 
to destroy a black hole via low energy proton bombardment. 

It is now an easy matter to determine the value of Q max that a given M can carry if it is to 
retain its character as a black hole. The maximum charge corresponds to that critical value, 
a c which makes V(A; M, Q max ) = E. For the case of E/2mc 2 = 1/400 (which corresponds to 
about 4.5 MeV in the case of a proton), we obtain a c ~ 9.6 x 10 -19 . Using this value in Eq. 
(TTSJ) yields the maximum value of charge that non-relativistic protons could add to a black hole 
of given mass, M: 

M. (19) 

IV. RELATIVISTIC PARTICLE DYNAMICS 

This section will extrapolate on the low-energy non-relativistic proton treatment. Consider 
now, depositing charge into a pseudo-Newtonian black hole by radial bombardment of high- 
energy protons. These particles are known to be abundant in astrophysical studies and to have 
a maximum observed energy of 10 20 eV.— Bombardment by N protons, each having charge q 
and mass m, increases the charge, Q, and mass, M, of the hole by N*q and N*m, respectively. 



Qmax ~ 8-3 X 10 



c_ 

kg 
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The charge addition to the hole significantly alters its character and interaction with additional 
bombarding protons, whereas the mass addition does not, and so it will be ignored. 

A particle having rest mass mo and velocity v while experiencing resultant force Fr will have 
equation of motion F R = ^(m 7f). The usual Lorentz factor, 7, is defined as 7 = ^_p 2 y/2 
and (3 = v/c. A proton radially approaching a charged pseudo- Newtonian black hole would 
experience competing gravitational and Coulomb forces, each of which is conservative, resulting 
in an equation of motion: 



d_ 

dr 



m O 70 + 



qQ 



(20) 



Using v = dr/dt, multiplying through by a differential dr and performing some elementary 
differentiation will show that: 



-d 



m O 70 + 



qQ 



moc 2 /3 
(1 - f3 2 ) 3 / 2 



dp. 



21) 



Amor 

This is easily integrated from infinity, the proton's presumed point of origination, to r, giving 

qQ 



m c [y r - 70]. 
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4:7TE r 

where 70 is the Lorentz factor for the proton at infinity and 7 r is its value at r. 

This can be solved for the Lorentz factor, 7 n for the proton at any r. 7 r will always have 
a value greater than or equal to 1. For a proton instantaneously at rest, such as at a turning 
point, Rt, a helpful inequality follows 



7o 



> 



1 j. 1 

— <Pr + - - 



qQ 



(23) 



This leads to a clear picture of the dynamics of a proton approaching the black hole. Suppose 
the proton has total energy at infinity, E = 7 m c 2 . Designate the right hand side of Eq. 
as F(r; Q, M) with given in Eqn Defining 

1 qQ 



F(r; Q, M) =2 In 



cosh I I)— 

r 



aM . Q 
—— smh — 
Q V r 



+ 



m c 2 Attest ' 

it is easy to find the maximum of F(r; Q, M) occurring at r = r ; the result is 

Q' 



tanh b 



r 



J^_p - 1 

r - _2_ ' 

L aM 



(24) 



(25) 



where T is a constant, defined as T = For protons T ~ 10 18 , and in order that the black 

amo 

hole exists, Q/aM < 1. Hence, to a highly accurate order of approximation, Eq. (1251) can be 
expressed as 



tanh ( b 



Q_ 

r 



Q_ 

aM 



1 

r' 



(26) 



A rescaling is convenient for graphing purposes. With again, A = GM / 2c2 ? an d a = ^ ■ This 
gives, from Eq. (l26il . 

tanh (y) = a " f • ( 2? ) 
The condition a < 1 in order for a black hole to exist, is demanded by the charge to mass ratio, 
Eq. f[T2"j) . Even in the case of a hole possessing a huge amount of charge, it will be the case 
that a <C 1. This allows for the approximation 

A min « — — f- (28) 

This equation represents a seemingly close approximation of distance of closest approach to the 
hole. The minimum A here is independent of the original speed given to the protons, because 
this equation represents the solution to the maximum of F(r; Q, M) for a certain charge-mass 
ratio, o. This maximum value of F(r; Q, M) correlates with the closest approach to the black 
hole for protons of a certain initial speed. 

Expressing F(r; Q, M) in its rescaled from, F(X, a), yields 



F(X,a) = 2 In 



cosh [ — ] sinh 

XJ a V A 



+ 21^. (29) 
A 



For highly relativistic protons, A will be approximately equal to 1, for the maximum value of 
F(X, a). This indicates that indeed the protons get very close to the black hole before turning 
around. Graphical analysis shows that A ~ 1 + , and implies that F(X, a) has a maximum turn 
around point when A > 1. Therefore the protons turn around before falling past the event 
horizon. To see this algebraically, that indeed A ~ 1, from Eq. ([9]) it can be shown that the 
black hole has a radius defined when the logarithmic argument goes to zero: 

tmh - lr (30) 

Using Eq. (1261) combined with Eq. (!30l) it can be seen that r is approximately equal to Rbh, 
but that ro is still slightly larger. These are related as, remembering that for protons, T ~ 10 18 , 

tanh ( bQ) = tanh ( b-^-] - -. (31) 
V r J V Rbh J T 

Now that we know the turning point will occur very close to the hole, for high energy protons, 

we can find the maximum value of F(X, a). Using A = 1 and some trigonometric substitutions, 

it is a straightforward matter to show from Eq. (T2"9"|) that 

F m ax(X, a) w 2 In f y ) + 2Ta. (32) 



To show that an upper limit exists, all that is needed is to prove that F max (\, <r) can be 
numerically greater than (70 - 1) where 70 = 10 11 and continue to have o < 1 as given by the 
restriction in Eq. (1121) . This value of 70, as mentioned before, belongs to extremely high-energy 
protons from an 'ultra-high' cosmic ray source. 

To ask what the maximum value of a is when F max is equal to 70 - 1 will be to ask what 
the charge-to-mass ratio of the hole is when it is maximally charged. The natural log term will 
be dominated by 2Ya max so that 

F m ax ~ 2r<J maa; . (33) 

From Eq. (|23|) . and safely ignoring the minus one because 70 is so large, we have 

70 ps 2Ta 

max • 

(34) 

Using 70 = 10 11 yields a max ~ 4.5 x 10~ 8 with T = q/arriQ for a proton having value ~ 1.1 x 10 18 . 
The value of a is also found graphically by plotting F(X, a) and 70, and finding their intersection. 
This means that all values of 70 for protons less than the energy given to them by 'ultra-high' 
cosmic rays will be stopped by a black hole with a max value of 4.5 x 10 -8 . Un-scaling, will give 
the charge-mass relationship 



Q « 3.9 x 10 



-18 



c_ 

kg 



M. (35) 



V. COMPARING LOW-ENERGY WITH HIGH-ENERGY 

It is helpful to see a how a comparison between low-energy protons and high-energy protons 
affects the upper-limit of charge. Low-energy protons with velocity 0.1c possess energy of about 
4.5 MeV, and an upper-limit of charge is given by 

Q « 8.3 x 10" 29 — M (36) 
_kg_ 

Compare this to Eq. (135]) and note that the order of magnitude of the constant is much 
less for the non-relativistic protons. It follows that with decreasing energy, the less capable the 
protons are of overcoming the Coulomb force. The less energetic protons will not be able to 
saturate the black hole as fully as the highly energetic protons. This explains the reason for a 
much smaller charge to mass ratio for low-energy protons. 

Note that the use of conservation of energy for the calculation of non-relativistic protons 
yields an equation of exactly the same form as the calculation using force for the relativistic 
protons. Equation (|24|) will be familiar in both calculations and set equal to E/mc 2 for the 
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low-energy protons, and equal to 70 — 1 for the high-energy case. Setting E/mc 2 = 70 — 1, it 
is easy to check that indeed they are the same equation. For a 70 value of 1.005 (with speed of 
0.1c) and using .5mv 2 for E (with again, the same speed of 0.1c) yields 0.005 as the answer for 
both. 

It can also be checked with high accuracy that Eq. ff28]) mirrors the low-energy calculation, 
which utilized conservation of energy. Using Ao ~ ojip — with a very precise o value of 
9.635104 x 10~ 19 yields 14.81462 for A . Using the low-energy calculation E/mc 2 = F(r; Q, M) 
with the same a value, yields 14.81433. The two results require a precise a, but indicate that 
both methods of determining Ao are accurate. The level of accuracy can be further increased 
with a more precise a value. 



VI. COMPARING PSEUDO-NEWTONIAN WITH GENERAL RELATIVITY 

Understanding how well pseudo-Newtonian gravity conforms to general relativity requires 
a comparison. Using the radial 'effective potential' equation for a charged particle falling in 
the Reissner-Nordstrom geometry^, we obtain the condition for maximal charge on the black 
hole when the bombardment of the last proton is unable to reach inside the event horizon. A 
turn-around point, where dr/dr = at r = r + , yields the condition, in raw form via natural 
units 

qQ = E(M+ ^M 2 - Q 2 ) . (37) 
Where E = '-fom, the energy of the proton at infinity. Solving for the charge to mass ratio, 

9- = 2E 1 ™ 

M q l + E 2 /q 2 ' 1 ' 

Converting units to SI, and making a simplification because Y = q/am w 10 18 3> 70 ~ 10 11 , 
yields 



Q 2 7o 



MVte^G r 

Rescaling gives 



(39) 



VGR = 2 -Y- (40) 

Comparing this result with the previous result, 70 ~ 2Ya pn) it follows that 

<?GR = 4*Tpn- ( 41 ) 
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It could be assumed that because the foundations of general relativity rely on its interpretation 
of physical geometry, and because space-time curvature was not taken into account for pseudo- 
Newtonian gravity, that the two descriptions in this case would differ drastically However, not 
only are they within the same magnitude, pseudo-Newtonian gravity is precisely four times too 
small. The curvature effects present in general relativity may completely explain the factor of 
four. Note especially that the pseudo-Newtonian radius for an uncharged black hole, derived by 
Peters^, GM/2c 2 , is also 4 times smaller than the Schwarzschild radius for an uncharged black 
hole, 2GM/c 2 . Because the charges on these black holes are far below the Reissner-Nordstrom 
condition, yet highly charged so that no known protons can get in, these radii are excellent 
approximations for the black holes involved. 

VII. CONCLUSIONS 

A merging of relativistic ideas with a non-relativistic theory yields an enlightening nonlinear 
theory of Newtonian gravity. An upper limit of charge was found for a black hole and it 
was made explicitly clear that a black hole cannot be destroyed by forcing protons into it 
as described by pseudo-Newtonian gravity. General relativistic predictions describe the same 
upper limit, and remarkably, the two results only differ by a factor of four. 
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